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The Tavis-Cummings model (the Dicke model treated in the rotating wave approximation) de- 
scribing many two-level systems coupled to a single bosonic mode, has been long known to show 
collective semiclassical oscillations when prepared in an inverted state, with all two-level systems 
excited, and the bosonic mode empty. This paper discusses how the quantum dynamics approaches 
this semiclassical result for large numbers of two-level systems, focussing on how the eigenvalues 
approach their semiclassical limit. The approach to the semiclassical result is found to be slow, 
scaling like a power of the logarithm of the system size. Considering also the effect of weak detuning 
between the two-level system and the bosonic field, quantum corrections are again found to decay 
slowly with system size, such that for a fixed detuning, the quantum effects of detuning are greater 
than the classical effect. 

PACS numbers: 42.50.Nn, 42.50.Pq, 03.75.Kk, 03.75.Lm 
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I. INTRODUCTION 

The Dicke model describing interaction between a 
number of two-level systems and a single bosonic mode, 
has long been studied as a simple model of cavity quan- 
tum electrodynamics, which despite its simplicity can 
show quite intricate behaviour. When the coupling be- 
tween the two-level system and the bosonic mode is 
treated in the rotating wave approximation, the Dicke 
model reduces to the the Tavis-Cummings model 0[H|- 
One long studied feature of this model is collective oscil- 
lations that arise when the initial state is inverted; the 
simplest such case concerns an initial state with an empty 
bosonic mode, and all two level systems in their excited 
state. In a semiclassical approximation [3[ the number 
of bosons describes a train of hyperbolic secant pulses. 
Interest in these collective oscillations has recently been 
revived both by the connection to atom-molecule inter- 
conversion in cold atomic gases 0, H, [E 5 as well as 
potential experiments studying coupling between quan- 
tum dots and a cavity photon modes [g, [9[ . The possi- 
bility of coupling between a radiation mode and multiple 
two-level systems is also being pursued in Circuit-QED 
experiments [10 ], in which the two- level systems are su- 
perconducting qubits. Both these latter examples are 
closely related to the original context of this problem Q : 
two-level atoms coupled to a photon mode in a cavity. 
This problem is also closely related to collective super- 
fluorescence (ll| for initially inverted atoms but without 
a cavity. Without the cavity there is only a single hy- 
berbolic secant pulse since the dense spectrum of photon 
modes prevents recurrence. 

The aim of this article is to study the quantum dy- 
namics of the Tavis-Cummings model with a finite num- 
ber of two level systems N, starting from a fully inverted 
state, in order to see how the quantum dynamics differ 



'Electronic address: jmjk2@cam.ac.uk 



from the semiclassical dynamics. In particular, consider- 
ing the case of a completely symmetric Tavis-Cummings 
model (where all two-level systems are identical), one 
finds that the semiclassical results are recovered in the 
limit of an infinite number of two-level systems, but that 
the approach to this semiclassical limit is slow, scaling 
as a power of the logarithm of the number of two-level 
systems. In addition to the sequence of hyberbolic se- 
cant pulses that exist in the semiclassical dynamics, the 
quantum dynamics is found to have an additional slow 
envelope. As the number of two-level systems increases, 
the period of this envelope increases compared to the 
period of the train of hyberbolic secant pulses and so its 
effects become negligible, however this trend is very slow, 
with Tcnvdopo/Tpuisc ~ [In ( y/N) ] 3 . 

A closely related question has been addressed by Farib- 
ault et al. [12], for the case of the Richardson model 
rather than the Tavis-Cummings model. Their work fo- 
cussed on how the integrability of the quantum model 
allows one to calculate overlaps between the initial state 
and the eigenstates, as well as matrix elements of the 
physical observables. The current work addresses a com- 
plementary question, that of how the eigenvalues ap- 
proach the semiclassical results for large system sizes. 

Aspects of the behaviour of the Dicke model in the ab- 
sence of the rotating wave approximation have also been 
studied; in particular, features of the finite size system as- 
sociated with the quantum phase transition in the infinite 
system size model have been considered. These include: 
changes of statistics of excited state energies [l3|, flij : 
perturbative approaches in the limit of large coupling 
strengths [TEl ]: entanglement between the bosonic mode 
and the two-level systems (l6| . including how that en- 
tanglement scales with the system size; and the scaling 
of other quantities, such as ground state energies or ex- 
citation gap with system size [l7], [HI- The results of 
these previous studies differ from the question addressed 
in the current work for two reasons: firstly, without the 
rotating wave approximation, the number of excitations 
in the system is no longer conserved, and so the classi- 
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cal problem is no longer integrable. Secondly, the results 
in this article relate to the collection of eigenstates near 
| -Eg | =0, while the ground state, or the thermodynamics 
at low temperatures (l9j , involves eigenstates with much 
lower energy. 

An alternate method to include quantum corrections 
to the semiclassical dynamics is by accounting for the 
dynamics of higher cumulants (as well as expectations) 
of the collective operators has been discussed by Vardi 
et al. Hl[ in related but different models. A simi- 
lar idea has also been discussed in the context the BCS 
model 0, HH, but in cases where the semiclassical dy- 
namics is more complicated. Another related problem 
concerns quantum dynamics in the central spin model — 
where a large nuclear spin (analogous to the bosonic field 
here) is coupled to a sea of electronic spins (analogous to 
two- level systems); Ref. 23] considers the quantum dy- 
namics of this model with initial conditions such that 
there is only a single excitation in the system. Other re- 
lated work concerns the dynamics of the Tavis-Cummings 
model in the opposite limit, of a small number of two- 
level systems, starting from an initial coherent boson 
state, studied in Refs. [13, HE HH in connection to the 
collapse and revival of Rabi oscillations in the case of a 
single two-level system. 

The results presented in this paper refer mainly to the 
completely symmetric Tavis-Cummings model, for which 
all two level systems are identical. This restriction al- 
lows one to extract simple analytical formulae for the 
scaling of quantum corrections with system size. In ad- 
dition, this symmetric case can be expected to be the 
"most classical" limit of the Tavis-Cummings model, as 
one may derive the semiclassical dynamics in such a case 
by replacing large quantum spins by classical spins. The 
fact that quantum corrections exist in this most classical 
case suggest that important corrections may exist in the 
non-symmetric Tavis-Cummings model. One specialised 
limit of this is discussed in Sec. IIV A[ supporting this 
idea. 

The rest of the paper is organised as follows; section llTl 
introduces the Hamiltonian, and discusses the previously 
known results of the semiclassical approximation; these 
are compared to the results of exact diagonalisation in 
Fig. [TJ Section IIIII then shows how the quantum cor- 
rections can be extracted from a WKB approach to the 
problem, focussing on the case where the two-level sys- 
tem and boson energies match; the effect of detuning is 
discussed in section lTVl and concluding remarks are given 
in section IVl 



II. THE MODEL AND COMPARISON OF 
SEMICLASSICAL AND NUMERICAL RESULTS 

The model studied in this paper can be written: 

N 



H = J2 M + 4 a + s i at ) . 



(1) 



where the spin operators obey [sf , sf] = ±sf , [sf, s~] = 
2s\ and a, are bosonic operators. The coupling be- 
tween the two-level systems and the bosonic mode has 
been scaled to 1, hence all other energies and times are 
measured in units of this coupling. The initial state of the 
system is taken to be \n = 0, ttt ■■■)■> where all the two- 
level systems are excited, and the bosonic mode is empty. 
In the subsequent dynamics there are collective oscilla- 
tions, transferring excitations between the two-level sys- 
tems and the bosonic mode. 

For comparison to the exact dynamics, the following 
briefly summarises the semiclassical solution, described 
in Refs. 0, The semiclassical equations correspond 
to writing the Heiscnbcrg equations of motion for the 
operators s~, sf, a, and then replacing these operators by 
commuting classical variables. The resultant equations 
can be solved by the ansatz: 



(ei — /i)a + ia 



(e; - /i) 2 + A 



(2) 



along with the equation of motion for a, d 2 = a 2 (A — a 2 ), 
and the self consistency conditions for /i, A: 



i 



(ei-/i) 2 + A' 



M = 



E 



(a -m) 2 + A' 



(3) 



These equations are analogous to the BCS gap equa- 
tion, with /i acting as a generalised chemical potential 
(i.e. common oscillation frequency), and A acting as the 
square of the gap (i.e. pairing field). The solution of the 
equation for a gives a train of hyperbolic secant pulses 
with a period T pu i so = 2 ln(\/A) / VA; such pulses can be 
seen in the time dependence of physical observables such 
as the occupation of the bosonic mode (n p hot} = a 2 - For 
the special case of ei = 0, Eq. © has the solution fi = 
and A = N. 

The problem this paper addresses can be seen most 
clearly in Fig. [T] which shows the dynamics of the pop- 
ulation of the bosonic mode n p hot according to Eq. (fT]) 
with ei = 0. The semiclassical train of hyperbolic pulses 
is seen in Fig. [TJ but there is in addition a slow enve- 
lope, not predicted by the semiclassical equations, and it 
is this slow envelope which is discussed below. 

By writing the photon number as a sum over eigen- 
states: 



i(t) = J2(o\x P )(x p \h\x q )(x q \o) t 



i(E p -E q )t 



, (4) 



one may note that the semiclassical result, with its per- 
fectly periodic train of pulses, corresponds (for e = 0) 
to: 



Eg = qSl, n 



N 



ln(ViV) 



(5) 



The inset of Fig. [T] shows the Fourier transform of 
n p hot{t)- In contrast to Eq. (O, the eigenvalues of the 
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full quantum problem are not equally spaced, and so 
the sidebands seen in Fig. QJb) arise. The period of the 
slow envelope is given by the splitting of these sidebands. 
Hence, to describe the approach to the semiclassical re- 
sult, one must consider how the deviation from regular 
spacing E q — q£l (i.e. the anharmonicity of the exact 
eigenvalues) collapses as N — > oo. 

Describing the quantum dynamics by finding the eigen- 
states is considerably simplified, since the initial state 
overlaps with only a small number of eigenstates. As 
discussed in Ref. [12], the initial condition chosen here 
satisfies (0|fl"|0) = 0, (0|# 2 |0) = N, while the semiclassi- 
cal eigenstates have E q = qir^/N / hi(y/N), which implies 
that only 0(ln(y/N) 2 ) eigenstates near \E q \ = can have 
any significant overlap with the initial state. The valid- 
ity of this is confirmed by noting that the result of exact 
diagonalisation and the result restricting the summation 
in Eq. Q to the seven smallest values of 
tinguishable by eye on the scale of Fig. [1] 



\E q \ are indis- 




FIG. 1: (Color online) Quantum dynamics of the number of 
photons for 2000 spins. Panel (a): time dependence of pho- 
ton amplitude. Panel (b) Principal Fourier components of the 
time dependence. The lowest seven eigenvalues are sufficient 
to describe the time dependence to better than visibie reso- 
fution on this scale. All energies and times are in units of the 
coupling between two-level systems and the bosonic mode. 



III. WKB APPROXIMATION AND SCALING 
OF CORRECTION 

To find the eigenvalues of the quantum problem, one 
may approach the problem by a method closely related to 
that of Bonifacio and Preparata [|[. As discussed above, 
this paper considers the symmetric case, = e, for which 
the Hamiltonian becomes H = eS z + S~a + S + a jf , where 
S = Si, and as a result, the quantum state may be 
written as a wavefunction in the one dimensional space of 
occupation, as used in Ref. Q to show how the semiclas- 
sical limit can arise. To find not only the semiclassical 
limit, but also the corrections to it, one may solve this 



one-dimensional problem using a discrete WKB approx- 
imation Using the basis \n) — \n phot = n,\S\ — 
N/2, S z = N/2 - n), the equation (E - eN/2)^ = 
becomes: 



(E + en)ip„ = nVN + 1 - nip n -i + {n+ N - mp n+1 

(6) 

The WKB approach consists of two parts; finding the 
WKB form of the wavefunction for (n, N — n) 3> 1, and 
then matching this wavefunction to appropriate forms 
for n ~ and n ~ N. In the following, this matching 
is referred to as matching the "boundary conditions" for 
the wavefunction at n = 0, n = N, but these boundary 
conditions are just the Schrodinger equation in Eq. ([B]), 
evaluated at n — 0, n = N, for which those approxi- 
mations valid for (n, N — n) 3> 1 no longer hold. This 
section will consider the case e = 0, the effect of non-zero 
e is discussed in section HVl 



The WKB form of the wavefunction can be derived by 
considering a transformation ip n — > i() n (i) n which gives 
the right hand side of Eq. ([6]) as a discrete derivative with 
a variable prefactor, i.e. Etp = —iv(x)d x ^p(x). Written 
this way suggests a WKB wavefunction of the form [27j 
tp w [1/ y / v(x)] exp[j J dzE/v(z)], or the discrete analog 
of this wavefunction. 

To write the wavefunction compactly, it is useful to 
introduce the function defined in Ref. [3( : 



9k 



V2r(l + A:/2) 




r[i/2 + fc/2] 



(7) 



y/k+ 1/2 fc»l 



which is chosen such that gkgk+i = k + 1. With this 
notation, the WKB wavefunction can be written: 



(WKB) 



cos(E$ n + (j) + nir/2) 
\Zgl9N-n 



(8) 



where § N = Y,Z=n 1 / 2 9? n 9N-m, and is an arbitrary 
phase set by boundary conditions. This wavefunction is 
valid while the change of phase between two successive 
values of n is much less than one; this condition for the 
validity of the WKB wavefunction can formally be writ- 
ten as 9n9N-n ^ E. Under the same conditions, one 
may approximate the summation in the definition of <£"„ 
by integration, and hence to leading order one has: 



f 



: cosh" 1 



»+3 



(9) 



It is worth noting that for even N, the solution in Eq. ([S]) 
with E = can be seen to be an exact eigenstate of 
Eq. @. 
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A. Matching WKB solution at boundaries 



can be written: 



1. Boundary condition at n — N 



When considering the boundary conditions, the bound- 
aries at n = and n = N behave diff eren tly. At 
n = N first note that g 2 N g ~ {N + l/2)y/2/n > E ~ 
iry/lV I ln(VA), meaning the WKB wavefunction is valid 
right up to this boundary. In addition, the Schrodinger 
equation at the boundary becomes EipN = NipN-i, 
and so to order 1 / \/]V, this boundary condition is sat- 
isfied by ipN-i = 0. Since the corrections to the semi- 
classical energies found below are of order l/\n(^/N), 
the much smaller corrections of order 1/yN can be ne- 
glected, and so the boundary condition at n — N requires 
4> — (2q + 2 — N)ir/2 where q is any integer; for even N 
one may thus choose <j) = 0. From herein, even TV is as- 
sumed; for odd N similar results with (f> = -k/2 follow 
straightforwardly. 



2. Boundary condition at n = 

The boundary condition at n = is more complicated. 
For this boundary, gfigN — (2/-k)^/N + 1/2, and so the 
condition g^gN E ~ ny/N /ln(-\//V) is not necessarily 
satisfied — the validity of the WKB wavefunction de- 
pends on taking ln(\/]V) ^> 1, which requires very large 
N; these corrections due to finite ln(viV) are the quan- 
tum corrections of interest in this paper. 

Because the WKB wavefunction breaks down for small 
values of n, it is necessary to match the WKB wavefunc- 
tion to the exact wavefunction at a non-zero value of 
n = no, rather than at n = 0; i.e the exact wavefunction 
should be calculated for n < no, and the WKB wavefunc- 
tion made to match it at n = no. The larger the value 
of no one takes, the better the calculated energies will 
match the exact solution. Surprisingly, even matching 
the solution at n — turns out to provide useful infor- 
mation, and correctly reproduces the scaling of energy 
with system size; this is discussed below in Sec. IIIIBI 

The procedure of matching is straightforward; for a 
given value no, one finds the exact solution ?Pn<n {E) by 
solving Eq. ©, with ip no +i = W™+\ B ^) i.e. 



"0 + 1 



v 2 \(x a \n )\ 2 , 

x g no 9N-n E _ x — • ( n ) 



This gives a nonlinear equation for E, the complexity 
of which increases with increasing no- An example of 
solving this equation, with no — 2 is shown in Fig.rjjb), 
which clearly accurately reproduces the results of exact 
diagonalisation. 




(27311 ln[2VN] 



FIG. 2: (Color online) Anharmonicity, E q — qE q , plotted 
against system size N with variables transformed to demon- 
strate agreement with the scaling expected from Eq. (|13|l . 
Panel (a) compares E q for q = 2, 3, 4. Panel (b) compares 
WKB approximations to E2 — 2E\ , on the same transformed 
axes, with no terms near the n — boundary. 



B. Simplified boundary condition and scaling with 
system size 

To understand the how the corrections to E q depend 
on number of two level systems N, the equation with 
general no is rather complicated, but as seen in Fig.[^b), 
matching at n = gives the correct scaling with system 
size, but an incorrect numerical coefficient. Considering 
this matching at n — 0, the boundary condition becomes 
from Eq. ©, ©: 



Ecos(E$ ) 



JVsin(E$i). 



(12) 



(10) 

where ff( n °) is the Hamiltonian restricted to n < no; the 
matching condition then becomes ipj$ = ip^ KB \ By 
diagonalising the (no + 1) x (no + 1) matrix to give 

ij("o) = A Q |x Q )(x Q |, the boundary condition equation 



In Fig. [2Ub), the solution of this equation is represented 
by the blue crosses, labelled no = 0. 

To the same level of approximation made so far, one 
may set $ - $1 = \n(VN)/\/N. If ln(\/F) were 
large, the solution of Eq. (fT2]l would be E q = qn/^i = 
qit\/~N / ln(V~N), i.e. the semiclassical result. For finite 
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ln(y/~N), Eq. (|T2|l produces corrections 0(q 3 ) that break 
the harmonicity; these terms look like: 



5E q = Cq 3 



N 



MVN)} 



(13) 



The coefficient C as calculated by expanding Eq. (|12jl in 
the small parameter l/hx(y/N) is C = tt 6 /24 ~ 40.1 
while the best fit coefficient over the range shown in 
Fig. (Ha) is C — 13.5, howevever this gradient may should 
be treated with caution, as 1/ ln(Vl0 6 ) is not a small 
number. The form in Eq. (|13p is the first main result 
of this paper — as shown in Fig. [2l this dependence on 
]n(yN) well describes the scaling of the anharmonicity; 
because the dependence on N is so slow, quantum cor- 
rections to the semiclassical result remain relevant even 
for 10 6 two-level systems. 




FIG. 3: (Color online) Time dependence of photon number 
for differing values of detuning e, showing modification to the 
slow envelope. 



IV. DETUNING 



Using the WKB wavefunction of Eq. (jSJ) one has: 



The remainder of this paper discusses the effect of de- 
tuning, i.e. of ej = e 7^ 0; this again reveals a distinction 
between the semiclassical dynamics and the full quantum 
mechanical problem, and provides some insight into the 
case where not all values of t\ are identical. The time 
dependence in the case of ti = e =/= is shown in Fig. [3j 
showing a change to the slow envelope. Semiclassically, 
the effect of non-zero e in Eq. ([3]) is to give /i = e/2, 
A = N — (e/2) 2 . Since it is A that controls the collec- 
tive oscillations, the leading correction to the frequency 
of pulses, i.e. the semiclassical energy splitting f2 should 
be quadratic in e: 



N 

71=0 



[l + (-l) n COS(2E,$ n )] 



9n9N-n 



N 

E- 

71=0 



l + (-l)"cos(2£ g $»)] 

9n9N-n 



(16) 



To analyse how the above expression depends on the 
number of two-level systems N , it is convenient to write it 
as (n q ) = (Aq + SA q )/(B + SB q ), where 8A q , SB q come 
from the term proportional to (—1)" in the summand, 
and A , B from the other term. One may then show 
that: 



E q = qSl, 



8VNln(VN) 



(14) 



The following discussion focusses on the case of small e, 
and shows that quantum corrections give shifts of ener- 
gies linear in e for small e, such that for any finite N, the 
leading correction to E q is quantum, not semiclassical. 

It is straightforward to see from Eq. ([6]) that for small 
e, the linear peturbative correction to eigenstate q would 
be: 



E„ 



E q -e(X q \h\X q ), 



(15) 



and so the aim of this section is to find how (n q ) = 
(Xq\n\X q ) depends on system size N and eigenstate label 
q. One should note that a q independent shift of E q has 
no effect on the time dependence of photon number. In 
the semiclassical picture, such a constant shift just cor- 
responds to a change of /i, hence our interest is in the 
quantity 8E q — SEq. 

Figure BJb) shows that the WKB eigenstates matched 
to the exact solution for n < n a — 3 reproduce the result 
of exact diagonalisation for N = 3000, giving confidence 
that the WKB wavefunction can be used to determine 



JV 

E 



9n9N-n 



dn—p 

o y/N-n 



= 2VN (17) 



and from Eq. ((9|), one sees that 

B « —^cosh^iVm) w —=\n(2V2N). 
s/N VN 



(18) 



Because Ao ^> Bo, one may anticipate that SB q has a 
more significant effect on the answer than does SA q ; this 
can be made more firm by noting that for even N: 



5A 



9=0 



(-1)" 



E 



4)" 



N 

y 



0.944, 



and numerically one may confirm that 5A q is a factor 
\/~N smaller than Aq in general. For 5B q one may observe 
that the main contribution to the sum comes from values 
of n <C N, where the denominator is smallest, and so 
approximate: 



SB a 



^E 

JV 71=0 



(-l)"cos(£ g $ r 
(n + 1/2) 



(19) 
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This expression for 8B q has the form of 1/v N multiplied 
by a function that should depend only on ]n(y/~N) — by 
considering the forms of E in Eq. and in Eq. ([5]), 
it is clear that only logarithmic dependence on system 
size enters the product E$ n . 

Putting all these considerations together one may 
write: 

2N 



2\n(VN) + 
and thus one may expand: 
2N 



NSB„ 



21n(VJV) 



NSB a 



ft 



ft 



ln(ViV) ln(VN) 



(20) 

for large N. This scaling is indeed seem to occur in 
Fig. (Ha) , and strongly suggests that /3° is independent 
of q, while f3q depends on q. One may thus write the q 
dependent part of the energy shift as: 



SE„ - 6E 



eN 



ft ~ ft 



ln{VN) V 21n(v7V) 2 



(21) 



where only the leading order term in powers of the log- 
arithm have been retained. This expression should be 
compared to the classical correction in Eq. (fH| . 




l/ln[2VN] 

FIG. 4: (Color online) Panel (a): Scaling of (n q ) with system 
size calculated using the WKB wavefunction matched to the 
exact solution for n < no = 40, and transformed following 
Eq. (|20|l . Dotted lines are a quadratic best fit to the data for 
N > 10 4 , assuming all lines tend to a common asymptote. 
Panel (b): comparison of exact diagonalisation and WKB ap- 
proximations for in q ) for N = 3000, showing that the WKB 
approximation should be valid for the panel (a). 



Comparing Eq. (f2Tj) to the semiclassical correction of 
Eq. (H3J, one may note that in both the limit fixed TV, e — > 
and also for fixed e, N — » oo, the quantum correction 
is larger than the semiclassical correction. 



A. Relating detuning to disordered t%. 

The above results on treating e perturbatively can also 
describe properties of the model where e;=i...jv-i = 0, 
epj = e; i.e. 5H — ts z N , again considering e perturba- 
tively. The fact this case can be treated by first or- 
der non-degenerate perturbation theory is not trivial; it 
arises because although degenerate states do exist (par- 
ticularly states with zero energy), the asymmetry intro- 
duced by altering a single spin energy does not mix 
these degenerate states. With a greater number of dif- 
ferent €i this simplifying condition fails, and degenerate 
perturbation theory is instead required. 

For the special form of SH above for which non- 
degenerate perturbation theory is relevant, one may con- 
sider the the expectation of perturbation Hamiltonian 
between the boson number states, for which one may 
show that 



(n\s*\n') = 8 n 



1 



n 

N 



(22) 



which reproduces the results of Sec. [IV] with e — ► e/N. 
Such a result hints that quantum corrections in the non- 
symmetric model may become more significant, since for 
weak disorder of energies e^, such linear corrections due 
to quantum corrections will win over the quadratic semi- 
classical effects of disorder. However, direct calculation 
of the quantum dynamics is challenging, as the size of the 
Hilbert space explored in such a model is exponential in 
N, and if quantum corrections still vanish as powers of 
ln(v)V), one requires very large systems to study the 
asymptotic behaviour at large N. 



V. CONCLUSIONS 

In conclusion, the quantum dynamics of the symmet- 
ric Tavis-Cummings model starting from an initially in- 
verted state describes a train of hyberbolic secant pulses 
(as in the semiclassical result), but with an additional 
slow envelope. The slow envelope arises due to the an- 
harmonicity of the eigenvalues of the quantum problem; 
this anharmonicity reduces as the system size increases, 
but only logarithmically with the number of two level 
systems. In the absence of detuning, i.e. for e, = 0, 
the eigenvalues take the form E q ~ qny/lV / ln(viV) + 
Cq 3 VN/[HVN)] 4 . 

With a small detuning e, the quantum problem has 
a perturbative correction linear in the size of detun- 
ing 8E (x e7V/[ln(v / ]V)] 3 , while the semiclassical result 
has only corrections quadratic in the detuning 8E oc 
t 2 /[Vn hi(-/N)]. As a result, for either finite detuning 
and N — > oo or finite TV and e — * 0, the classical ef- 
fects vanish faster than quantum corrections. This result 
may also be indicative of the effects of disorder, i.e. of 
non-identical values of q. For a fixed distribution of e,, 
as N — > oo the effects of this disorder vanish compared 



7 



to the energy scale of the common coupling ~ >/~N. It 
seems likely that in such a limit, quantum corrections 
may vanish more slowly than the classical effects of the 
distribution of e^. 

In the general case of disordered ej, the quantum prob- 
lem is significantly harder to solve, as the size of the 
Hilbert space grows exponentially with the number of two 
level systems, whereas for = e, it grows only linearly. 
However, the integrability of the problem may simplify 
the problem, as discussed by Faribault et al. [13]. In 
the symmetric case, the integrability of the semiclassical 
problem means that while conservation laws restrict the 
system to exploring an iV dimensional Hilbert space, in 
the limit of large N, the system in fact only explores a 
one dimensional path through this space. Similarly for 
the disordered model, the problem is integrable [H, [H|, 



which in turn leads to the semiclassical dynamics [J, [H, 0] 
following a one dimensional path, suggesting that a sim- 
ple description of the quantum corrections even for large 
values of N might be possible. 

Note added: Since the submission of this work, a sim- 
ilar treatment of this model has been undertaken by Ba- 
bclon et al. [30(. 
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